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ABSTRACT 

I present a method to deproject the observed intensity profile of an axisymmetric bulge 
with arbitrary flattening to derive the 3D luminosity density profile and to calculate 
the contribution of the bulge to the rotation curve. I show the rotation curves for 
a family of fiducial bulges with Sersic surface brightness profiles and with various 
concentrations and intrinsic axis ratios. Both parameters have a profound impact on 
the shape of the rotation curve. In particular, I show how the peak rotation velocity, 
as well as the radius where it is reached, depend on both parameters. 

I also discuss the implications of the flattening of a bulge for the decomposition 
of a rotation curve and use the case of NGC 5533 to show the errors that result 
from neglecting it. For NGC 5533, neglecting the flattening of the bulge leads to an 
overestimate of its mass-to-light ratio by approximately 30% and an underestimate of 
the contributions from the stellar disc and dark matter halo in the regions outside the 
bulge-dominated area. 
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1 INTRODUCTION 

Rotation curves of spiral galaxies are one of the prime 
methods to probe the mass distri bution on galactic scales . 
The discovery, made in the 1970s dRogstad fc Shostakll 19721 : 
iRoberts fc Whitehursd 1 19751 : iBosmal Il978l : ; IBosmal ll98lF , 
that rotation curves remain flat till the last measured points, 
well outside the optical discs, provided the decisive evidence 
for the presence of large amounts of dark matter in galax- 
ies. Unfortunately, however, it has proven difficult to gain 
detailed knowledge about the distribution of dark matter in 
galaxies. This is primarily due to the fact that, with the ex- 
ception of the lowest surface brightness galaxies, the stars 
in galaxies contribute a significant, but not prec isely known, 
fract i on to the overall gravitational potential (|Persic et alj 
19961: Ide Blok fc McGaughll 19971 ; IPahmas fc Williamsll2000l ; 



luminosity-profile (|de Vaucouleursl Il958l : Ivan Houtenl Il96ll : 
Ide Vaucouleurslll974 V For dynamical purposes, they were 
often further assumed to be spherically symmetric bod- 
ies, in which case the observed intensity profile can be de- 
proje cted and the gravitational potential calculated easily 
(e.g. lYound Il976t ). In reality, however, bulges are gener- 
ally nei ther spherical, nor do they follow the R 1 ^ 4 density- 
profile l | Andredakis fc Sanders 1994; Andrcdakis et alJl995| : 
Erwin fc Sparke 1999; Graham 200li; lErwin fc SpaikelEiOOa 



Noordermeer fc van der Hulsdl2007l ). 
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To determine the contribution from the stars to an 
observed rotation curve, the full 3D structure of the stel- 
lar mass distribution needs to be known. In most disc 
galaxies, the bulk of the stars can be divided between 
a flat disc and a thicker central bulge. Once an accu- 
rate bulge-disc decomposition has been performed, it is 
relatively straightforward to calculate the contribution to 
the rotation curve from t he disc, given it s simple, effec- 
tively 2D structure (e.g. ICasertanol |l983). For the cen- 
tral bulge, however, the situation is more complex. His- 
torically, it was often assumed that bulges follow an 7? 1//4 
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Several papers have addressed various aspects related 
to the properties of non-spherical bulges with a more gen- 
eral luminosity density profile. The problem of deproject- 
ing a non-spherical bulge was already addressed by IStarkl 
l|l977h . The dynamical properties of spherical bodies follow- 
ing a g enera l 7? 1//n luminosity profile have be en studied by 
ICiottil (|l99ll ). This analysis was extended bv lTruiillo et al.l 
(2002), who studied triaxial systems following the R 1 ^ lu- 
minosity law and presented an analytical expression to ap- 
proximate the deprojected, 3D luminosity density profile. 

For the purpose of rotation curve analysis, one usually 
assumes axisymmetry for all components in the galaxy and 
ignores the triaxi al nature of th e bulg e. Most bulges are only 
mildl y triaxial (jBertola et al.l Il99ll : iMendez-Abreu et al.l 
l2007h . so that the errors stemming from this assumption 
are generally small. Accounting for the full triaxiality of the 
bulges merely complicates the calculations in this context. 
In fact, given the fact that rotation curves are by definition 
one-dimensional quantities, they are intrinsically unsuitable 
to account for deviations from axisymmetry. A proper treat- 
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ment of the dynamical effects of triaxial bulges requires two- 
or three-dimensional modelling of the gas flows in their po- 
tentials. The flattening of the bulge in the vert ical direction, 
on the other hand, is often significant (e.g. iBertola et al.l 
ll99ll ; lNoordermeer fc van der Hulstll2007t) . It is also dynam- 
ically important and should not be neglected, even in the 
axisymmetric approximation. 

In this paper, I study the 3D luminosity density of ax- 
isymmetric bulges with various intrinsic flattenings and their 
contribution to the rotation curve of a galaxy. In section [2] I 
present the relevant equations to deproject the observed sur- 
face brightness profile of a bulge to the 3D luminosity den- 
sity and to calculate its contribution to the rotation curve. 
Although some of the equations in section [2] can be applied 
to an arbitrary bulge surface brightness profile, I specificall y 
focus on bulges with a Sersic luminosity profile dSersidl 1968T ) . 
In section [3] I present the bulge rotation curves for a range 
of parameters and discuss various properties of the curves. 
I show the effects of the concentration and flattening of the 
bulge on the shape of its rotation curve and in particular on 
the location and height of the peak. Finally, in section [4] I 
briefly discuss the implications of the results from the pre- 
ceding sections. In particular, I use the example of the Sab 
galaxy NGC 5533 to illustrate the errors which result from 
neglecting the flattening of the bulge when decomposing a 
rotation curve. 




Figure 1. Schematic overview of the two coordinate systems. 
Black lines give the galactic coordinate system, grey lines give 
the coordinate system of the observations (see text). The ellipse 
shows an (x, z)-cut through an isodensity surface. 



2 DERIVATION OF THE ROTATION CURVES 

2.1 Deprojection of the observed surface density 
distribution 

The first step in our calculations involves the deprojection 
of the observed surface density distribution, to yield the 3D 
luminosity densit y. This has b een done for the general tri- 
axial situation by IStarkl(|l977r ). Here, I consider the simpler 
configuration of an axisymmetric spheroid, symmetric with 
respect to the plane of the galactic disc. 

Let us define an (x, y, z)-coordinate system such that 
the (x, j/)-plane coincides with the plane of symmetry of 
the galaxy. The luminosity density can then be written as 
p(x,y,z) — p(m), with m = ^/x 2 + y 2 + (z/q) 2 and q the 
intrinsic axis ratio of the bulge isodensity surfaces. 

A second coordinate system, (C, K ,£)i can be defined 
corresponding to the observational configuration. £ lies along 
the line of sight and is inclined by an angle i with respect 
to the z-axis. The re-axis coincides with the j/-axis; both 
are defined to lie along the line of nodes. £ is the axis that 
is perpendicular to re on the plane of the sky; in practice 
this means that £ lies along the minor axis of the galaxy. 
See figure [T] for an overview of the coordinate systems. The 
(x,y,z) and coordinate systems are related by the 

following transformation matrices: 
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position (ft, £) on the sky can be written, in the case of an 
optically thin medium, as 
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It can be shown that the projected intensity distribution on 
the sky has an ellipsoidal shape, where the ellipticity e of 
the projected image is related to the intrinsic axis ratio q 
and the inclination angle i as 



(1 — e) 2 = q 2 + (1 — q 2 ) cos 2 % 



(4) 



Thus, the problem is reduced to finding the relation between 
the projected intensity at the line of nodes (£ = 0) and 
the density p(m). Once the intensity at the line of nodes is 
derived, the rest of the image can be reconstructed using 
equation [4] 

Given a position fto at the line of nodes (f = 0), the 
relation between £ and m is (see equations [T] [2j : 



m 2 (C| fto,£ = 0) 



x 2 «| ft ,£ = 0)+y 2 «| ft ,£ = 0) + 
-V(C| ft f „£ = 0) (5) 



= £ sin i + ft + 



yi 2 • 
Q cos i 



C 2 (m| ft O ,£ = 0) 

K 

dm 



2 2 

m — k 



sin 2 i + 4r cos 2 i 



sin 2 i + cos 2 i ) (m 2 — re 2 ,) 



(6) 

m 



With these coordinates, the projected intensity at each 



Then, setting £ = in equation [3] changing the integration 
variable from £ to m and using equation [7] the projected 
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intensity I(ko) at the line of nodes becomes: 

/(«o) = / p(c|K ,e = o)dc 



= 2 



/ p(m)—dm 
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(8) 
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m dm 



This is an Abel-integra l, which can be inverted (e.g. 
iBinnev fc Trem ainc 1987, app. 1.B4) to obtain the following 
equation, which gives the 3D luminosity density distribution 
p(m) when the observed major-axis intensity profile I(k) is 
given: 



p(m) = ./sin i+ — cos 2 i / ___=_. (9) 

f V 9 J m dn^ K 2 -m 2 



2.2 Calculation of the rotation curve 

Now, assuming that the observed emission traces mass 
and that the mass-to-light ratio is constant throughout the 
bulge, equation|9]also gives the 3D mass density distribution 
in the bulge (except for a constant factor for the actual value 
of the mass-to-light ratio, T). The derivation of the corre- 
sponding rotation curve is then stra ightforward and is de- 
scribed in IBinnev fc Tremainel (1 19871 ). Inserting equation [9] 
into their equation 2-91, we get for the rotation curve: 



V 2 (r) = —AGqT 



' i H — - cos 2 
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where e is the eccentricity of the bulge: e = — q 2 . 

This equation is valid for any given observed intensity 
profile /(ft); no prior assumptions have been made regard- 
ing the shape of the light or density profiles, except that of 
spheroidal symmetry. In principle, one could directly derive 
the rotation curve for any given intensity profile by numer- 
ically calculating the derivative dl/dn and evaluating the 
integral. In practice, the inte nsity profiles of many bulges 
follow a general Sersic profile (Sersic 1968|), defined as: 



I(,(k) = loexp 



i/t. 



(11) 



Here, 7o is the central surface brightness and n is a concen- 
tration parameter that describes the curvature of the profile 
in a radius-magnitude plot. For n = 4, equation II lj reduces 
to th e well-known de Vaucouleurs profile l)de Vaucouleursl 
1 19481 ). whereas n = 1 describes a simple exponential pro- 
file, ro is the characteristic radius, which is related to the 
effective radius (r e , the radius which encompasses 50% of 
the light) as r e = b" ro. b n is a scaling constant that is de- 
fined such that it satisfies 7(&n, 2n) = ^r(2n), with 7 and T 
the incomplete and complete gamma functions respectively. 
Given a Sersic bulge, one can evaluate the derivative dl/dn 



analytically: 

— - 
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Inserting this into equation 1101 we get the following, 
final equation for the bulge rotation curve: 
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3 PROPERTIES OF THE ROTATION CURVES 

Equation [T3] can, to my knowledge, generally not be solved 
analytically. In stead, it was evalua ted numerically, using al- 
gorithms from Press et all l|l992lfl for various parameter 
combinations. 

Of particular interest are the effects of the flattening, 
q, and the concentration, n, of the bulges on the result- 
ing rotation curves. In figure [2] I show the rotation curves 
for a family of fiducial bulges with various concentrations 
and axis ratios, but with 7o and ro scaled such that all 
bulges have equal total luminosity (Lt = 5 • 10 9 Lq) and 
effective radius (r e = lkpc). For simplicity, I assumed a 
mass-to-light ratio of T = 1. The range of values for n and q 
which I consider here spans the observed values in bulges 
of nearby spiral galaxi e s of various morphological types 
llAndredakis et aL 19951: Khosroshahi et all l2000l : iGrahanj 
120011 : iNoordermeer fc van der Hulstll2007l) . although I also 
include a number of prolate bulges with q > 1 which are 
not observed in reality. Note that observed bulges show a 
general trend of increasing luminosity and radius with in- 
creasing n. For easy comparison, however, I choose to study 
bulges with equal luminosity and effective radii here; the 
rotation curves shown in figure [5] can be scaled easily for 
bulges with different mass or radius. 

It is clear that both the flattening and the concentra- 
tion of the bulge have a large influence on the shape of the 
rotation curves. Varying either parameter changes both the 
radius where the rotation curve peaks and the correspond- 
ing peak velocity. This is illustrated in a quantitative way 
in figure |31 As expected, a flattened bulge has a higher peak 
rotational velocity than a spherical one. Interestingly, the 
strength of the correlation is a weak function of the con- 
centration, being largest for less concentrated bulges. For 
n — 0.5, a highly flattened bulge with an intrinsic axis ratio 
of q — 0.25 has a peak velocity about 22% higher than the 
spherical case. For a highly concentrated bulge (n = 8), the 
difference is approximately 17%. 

Not surprising either, more centrally concentrated 
bulges reach the peak in their rotation curve at smaller radii. 



1 A FORTRAN program to evaluate the integrals in equation [13] 
is available upon request from the author. 
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Figure 2. Rotation curves for bulges with total mass M;, = 5 ■ 10 9 Mq and effective radius r e = 1 kpc. The different panels are for 
different values of the Sersic concentration parameter n, as indicated in the top-right corners. The solid lines show, from top to bottom, 
rotation curves for axis ratios q = 0.2, 0.4, . . . , 1.0 respectively. The dashed lines are for the unrealistic cases of prolate bulges (q = 1.5 
and 2.0). The dotted lines show the Keplerian rotation curve for a point-like body of equal mass. 
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Figure 3. The dependence of the peak rotation velocity (Vpeak, top panels) and the corresponding radius (iJ pea ki bottom panels) on 
the bulge flattening (q, left hand panels) and the concentration parameter (n, right hand panels). Different colours and symbols in the 
left hand panels indicate different values for n, as labeled; dashed lines indicate unrealistic values of q corresponding to prolate bulges. 
The different lines in the right hand panels show the behaviour for different values of q. 
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Cored bulges with n = 0.5 reach the maximum rotation ve- 
locity between 1.5 and 2.0 effective radii, depending on q. 
Concentrated bulges with n = 8, on the other hand, peak at 
very small radii, within 2% of the effective radius, with the 
exact radius again depending on q. 

The dependence of the location of the peak in the ro- 
tation curve on the flattening of the bulge is somewhat less 
expected. The bottom left hand panel in figure [3] shows that 
the peak is reached at smaller radii in more flattened bulges. 
The strength of the correlation depends on the concentra- 
tion of the bulge: for n = 0.5, a highly flattened bulge with 
an intrinsic axis ratio of q — 0.25 reaches the peak rotation 
velocity at a 12% smaller radius than the spherical case, 
whereas for n = 8, this difference is 18%. As an aside, I note 
that a spherical bulge (q — 1) with a Sersic parameter of 
n — 2 reaches the peak rotational velocity at exactly one 
effective radius. 

Finally, for a given total mass, effective radius and in- 
trinsic axis ratio, the peak rotational velocity is lowest for a 
bulge with n = 1.5, with both more and less concentrated 
bulges having higher peak rotation velocities. 



4 DISCUSSION 

It is now common knowledge that morphologically, the 
assumptions that bulges are spherical and have an i? 1//4 
luminosity-profile are poor descriptions of reality. Our anal- 
ysis shows that this is also true dynamically. As figures[2]and 
[3] illustrate, the flattening and concentration of a bulge have 
a strong influence on the shape of its rotation curve. Specif- 
ically, non-spherical bulges and bulges with a concentration 
parameter n different from 4 have rotation curves with very 
different peak velocities and radii than a 'classical' spherical 
de Vaucouleurs bulge. 

The dependence of the peak velocity on the flattening 
of the bulge has important consequences in the decomposi- 
tion of an observed rotation curve into contributions from 
the various components in a galaxy. In practice, the bulge 
contribution in such a decomposition is often determined 
by scaling up the bulge mass-to-light ratio until the peak 
of the bulge rotation curve matches the amplitude of the 
observed curve at that radius. If we neglect the intrinsic 
flattening of a bulge and instead derive its rotation curve 
under the assumption that it is spherically symmetric, we 
will underestimate the peak velocity (see figure [2}. As a re- 
sult, we will need to assume a higher mass-to-light ratio in 
order to explain the observed rotation velocity at the peak 
radius. Figure [3] shows that the difference in peak velocity 
between a spherical bulge and a highly flattened bulge with 
q — 0.25 is of the order of 20%. Since the mass-to-light ratio 
is proportional to the square of the bulge rotation velocity 
(equation I13|) . the error in the former would be more than 
40%. 

To illustrate this effect, I show in figure [4] the ro- 
tation curve decomposition for the massive Sab galaxy 
NGC 5533. Its rotation curve was measured from a com- 
bination of_0£ticaI_s_p_ectoos£0£ic data and Hi observa- 
tions by iNoordermeer et al.l ( 2007 ). R-band photometry for 
this g alaxy was presented by Noordermeer fc van der Hulstl 
(2007), as well as a detailed study of the properties of its 
bulge and disc. The effects of the seeing on the measured 



Table 1. Results from the rotation curve decomposition of 
NGC 5533. The middle column shows the fitted parameters for 
the case where the flattening of the bulge was taken into account, 
while the right-hand column is for the situation where the bulge 
was assumed to be spherically symmetric. The errors are the for- 
mal fitting uncertainties, and do not include any systematic un- 
certainties, nor the degeneracies between the various parameters. 





q = 0.33 


q = l 


T 6 (M /L 0jJl ) 


2.8±0.1 


3.7±0.1 


T d (M Q /L 0iH ) 


5.0±0.2 


4.4±0.2 


Rc (kpc) 


1.4±0.2 


1.7±0.2 


po (M Q pc- 3 ) 


0.31±0.08 


0.23±0.05 


M BH (1O 9 M ) 


2.7±0.6 


3.4±0.6 



bulge par ameters were c orrected for, using the deconvolu- 
tions from lCrahaml (|200ll ) . The bulge was found to be highly 
flattened, q = 0.3, making this galaxy an interesting test- 
case for the effect described above. 

Based on the derived p r operti es of the bulge from 
INoordermeer fc van der Hulstl (|2007h . I have derived two 
versions of its contribution to the rotation curve: one ac- 
counting correctly for its flattening using equation [13] and 
another by neglecting the flattening and assuming q = 1. 
The contribution of the stellar an d gaseous d i scs we re calcu- 
lated using the prescriptions from ICasert and l)l983h . assum- 
ing a vertical scale height of one fifth of the stellar disc scale 
length. For simplicity, the dark matter halo was assumed 
to have an isothermal density profile. Finally, in order to 
explain the rapid rise of the rotation curve in the centre, 
a central black hole was added to the model. The contri- 
butions from the various components were then added and 
the 5 free model parameters (the R-band mass-to-light ra- 
tios T;, and Td for the bulge and disc respectively, the core 
radius R c and central density po of the dark halo and the 
mass of the black hole Mbh) were adapted in a simple \ 2 
minimisation procedure to obtain the best fitting model ro- 
tation curve. The values for the fitted parameters are listed 
in table [T] 

Figure [4] shows that equally good fits can be made with 
both bulge rotation curves; neglecting the bulge flattening 
does not have a noticeable effect on the quality of the model 
in this case. However, it is clear from table[T]that, when the 
bulge is assumed to be spherically symmetric, a higher mass- 
to-light ratio has to be assumed to explain the observed 
rotation velocities around the peak radius of 2 kpc. In this 
case, the difference between the two values for Tj, is 32%. 

Figure [4] and table [T] also illustrate a more subtle ef- 
fect. The effect of the flattening of the bulge on its rotation 
curve contribution diminishes rapidly with radius (see also 
figure [2]). Outside about 5 effective radii, the flattening has 
virtually no effect anymore on the amplitude of the bulge ro- 
tation curve, which then only depends on the total mass of 
the bulge. Thus, because we were forced to assume a higher 
mass-to-light ratio for the spherical bulge, its contribution 
at large radii is larger than in the flattened case. As a re- 
sult, the derived contribution of the disc is somewhat smaller 
and the dark matter halo is less concentrated, with a 21% 
larger core radius. The s ame effect occurs wh en we assume 
an NFW density profile (|Navarro et al.lll997T ) for the dark 
matter halo. In that case, we find again that the concen- 
tration of the halo is significantly smaller (by about 15%) 
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Figure 4. Rotation curve decomposition for NGC 5533. Top panels show the decomposition using the correct value of the bulge flattening. 
The bottom panels show the case where the bulge flattening was neglected and its contribution was calculated under the assumption of 
spherical symmetry. Left hand panels show the central regions while the right hand panel s show the full rotation cu rve. The data points, 
errorbars and gray shaded band show the observed rotation curve and uncertainties (see Noorder meer et alj|2007l for details). The thin 
lines show the contributions from bulge (solid), stellar disc (long dashed), gaseous disc (dotted), central black hole (short dashed) and 
dark matter halo (dot-dashed). The bold line shows the rotation curve for the combined model. 



when we assume that the bulge is spherical, compared to 
the model where we account for the flattening of the bulge. 

In conclusion, the results shown in this paper prove 
that, for a detailed study of the dynamical structure of disc 
galaxies with a sizeable bulge, it is crucial to properly ac- 
count for the concentration and intrinsic flattening of the 
latter. 
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